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1.0  SUMMARY 


Development  of  advanced  propulsion  systems  of  interest  to  the  Air  Force  requires  study  of  the 
harsh  p  hysical  environments  as  sociated  w  ith  reacting  f  lows  eh  aracterized  b  y  extremes  i  n 
temperature  a  nd  pr  essure  a  nd  c  omposition  i  n  t  he  c  ase  of  pi  asmas.  S  uch  e  nvironments  a  re 
difficult  to  s  tudy  in  d  etail  w  ith  me  thods  th  at  ma  ke  u  se  o  f  p  hysical  i  ntrusion  in  to  th  e  lo  cal 
environment.  Such  difficulties  motivate  the  use  of  non -contact  means  of  study  such  as  optical 
approaches. 


By  taking  advantage  of  impressive  advances  in  diode  laser  and  associated  technology,  tunable 
diode  1  aser  a  bsorption  s  pectroscopy  (  TDLAS)  ha  s  be  come  a  w  idespread  t  ool  f  or  opt  ical 
diagnostics  i  n  r  eacting  f  lows.  M  easurements  of  t  arget  s  pecies  c  oncentration,  t  emperature, 
pressure,  a  nd  ve  locity  h  ave  be  en  de  monstrated  i  n  a  w  ide  va  riety  of  e  nvironments  b  y  s  everal 
research  groups.  C  ommercialization  of  ni  che  m  easurement  i  nstruments  ba  sed  on  T  DLAS  h  as 
been  achieved.  The  current  maturity  of  the  technique  enables  its  application  to  hypersonic  flow 
paths  -  scramjet  e  ngines  i  n  pa  rticular.  In  t  his  r  eport  w  e  s  ummarize  our  pr  ogress  i  n  a  pplying 
TDLAS  t  o  s  cramjet  i  solators  a  nd  c  ombustors.  W  e  a  iso  p  erform  close  s  crutiny  o  f  c  ollected 
signals  to  ascertain  expectations  of  accuracy  and  precision  in  measurements  coupled  to  reacting 
flows.  W  e  ex  amine  s  ome  T  DLAS  d  ata  s  ets  i  n  th  is  c  ontext  a  nd  discuss  t  he  pr  opagation  of 
uncertainties  from  primary  measured  values  of  signal  amplitude  versus  time  to  secondary  values 
of  temperature,  number  density,  and  pressure. 


2.0  INTRODUCTION 


Application  of  tunable  diode  laser  absorption  spectroscopy  (TDLAS)  to  reacting  flows  has  been 
repeatedly  d  emonstrated  i  n  va  rious  e  nvironments  pr  imarily  fort  emperature  de  tennination. 1  At 
least  one  s  tudy  de  monstrated  s  ystem  c  ontrol  us  ing  T  DLAS-based  m  easurement  i  nput  i  n  r  eal 
time.  In  this  work,  closed-loop  temperature-based  control  of  a  pulsed  incinerator  was  achieved 
by  m  aximizing  th  e  me  asured  time  -averaged  t  emperature  t  hrough  adjustment  of  t  he  pha  se 
between  primary  and  secondary  air  forcing.  Overall,  the  TDLAS  technique  has  matured  to  the 
point  that  several  c ompanies  including  Zolo  Technologies  Inc.,  Southwest  Sciences,  Inc.,  Los 
Gatos  R  esearch,  Inc.,  a  nd  P  hysical  S  ciences  Inc.  h  ave  d  elivered  T  DLAS-based  co  mmercial 
products  for  niche  measurement  applications.  In-stream  measurements  of  temperature,  pressure, 
species  c  oncentration  a  nd  ve  locity  a  re  pos  sible  i  n  h  ypersonic  r  eacting  flows  us  ing  T  DLAS.1 
Specifically,  d  emonstration  m  easurements  oft  emperature  and  w  ater  co  ncentration  h  ave  b  een 
made  at  the  exit  of  a  d  irect-connect  model  scramjet  engine,  and  air  velocity  has  been  measured 
in  a  d  irect-connect  isolator  using  water  and  oxygen  absorption.4  Measurements  of  temperature 
and  w  ater  concentration  ha  ve  a  Iso  b  een  r  eported  i  n  w  ind-tunnel  te  sting  o  f  mo  del  s  cramjet 
engines.3  More  recently,  spatial  information  about  the  shock  train  in  a  scramjet  isolator  has  been 
captured  u  sing  mu  ltiple  T  DLAS  lin  es  o  f  sight  a  nd  a  c  luster  m  ath  a  lgorithm  a  pproach  t  o 
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tomography.6  These  results  are  all  highly  encouraging  and  point  toward  successful  integration  of 
the  t  echnique  i  nto  ground  a  nd  flight  t  ests.  T  he  H  ypersonic  International  F  light  R  esearch 
Experimentation  ( HIFiRE)  pr  ogram  ha  s  e  mbraced  T  DLAS  t  echnology.  The  F  light  1  pa  yload 
includes  t  wo  ox  ygen-based  i  nstruments  f  or  m  ass  c  apture  d  etennination.  A  nd,  w  ater-based 
measurements  of  concentration,  temperature,  and  pressure,  will  be  executed  on  Flight  2  using  8 
lines  of  sight  at  the  exit  plane  of  the  combustor. 


A  mu  ltiple  1  ine-of-sight  t  unable  di  ode  1  aser  a  bsorption  s  ystem  f  or  m  easuring  s  hock  t  rain 
structure  inside  a  model  scramjet  isolator  at  the  Wright-Patterson  AFB  Propulsion  Directorate’s 
research  cell  1 8  supersonic  combustion  facility  is  described.  Multiple  rovibrational  lines  of  water 
are  p  robed  us  ing  t  ime-multiplexing  o  f  mu  ltiple  d  iode  la  sers  a  long  1  6  lines-of-sight  th  at  a  re 
projected  across  a  rectangular  isolator  duct.  S  tatic  temperatures  and  water  concentration  (path- 
averaged)  a  re  c  omputed  f  rom  a  B  oltzmann  pi  ot  f  or  t  he  m  ultiple  1  ines  t  hat  a  re  m  easured,  a  nd 
static  p  ressures  a  re  d  etermined  from  t  he  a  bsorption  1  inewidth.  The  ex  perimental  s  etup  an  d 
measurement  technique  are  presented  here  along  with  the  Mach  2  flow  (no  shock  train)  results. 
The  absorption  measurements  agree  to  within  2%  of  sidewall  pressure  transducer  measurements 
and  calculated  (path- averaged)  static  temperatures. 


3.0  METHODS 


3.1  Measurement  Technique 


Because  of  the  ubiquity  of  the  use  of  absorption  spectroscopy  in  multiple  fields,  many  different 
notations  a  re  us  ed  t  o  describe  t  he  s  ame  ph  enomenon. 7,8,9  For  cl  arification,  t  he  p  articular 
formulation  us  ed  he  re  i  s  de  scribed  br  iefly  be  low  a  nd  f  ollows  of  t  hat  de  fined  i  n  t  he  H ITRAN 
database.10,11 

Beer’s  Law  describes  the  attenuation  of  light  as  it  passes  through  an  absorbing  medium:7 
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In  Equation  1,  /  is  the  transmitted  laser  intensity  through  the  medium,  Io  is  the  initial  intensity,  oy 
is  the  absorption  cross  section  at  the  frequency  v,  AN  is  the  difference  between  upper  and  lower 
state  q  uantum  s  tate  d  ensities  ta  king  in  to  a  ccount  d  egeneracy  factors  (  total  q  uantum  s  tate 
fractional  population  times  the  total  species  density),  L  is  the  path  length  through  the  absorbing 
medium,  and  av  is  the  absorption  coefficient  at  the  frequency  v.  The  absorbance,  A  „,  shown  in 
equation  1  i  s  pi  otted  i  n  this  report  rather  t  han  t  he  a  ctual  i  ntensity.  T  his  i  s  done  b  y  di  viding 
through  with  a  reference  signal  that  is  recorded  in  addition  to  the  measurement  paths.  However, 
Equation  1  i  s  not  i  n  a  convenient  f  orm  f  or  comparison  w  ith  w  hat  i  s  typically  t  abulated  i  n 
spectroscopy  databases  such  as  HITRAN.1 1  These  databases  replace  the  absorption  cross  section 
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with  the  lin  e  in  tensity,  S(T),  and  a  1  ineshape  f  unction  g(v-vio)  that  describes  t  he  s  pectral 
lineshape  w  here  vio  is  the  pe  ak  c  enter.  In  a  ddition,  t  he  1  ine  i  ntensity  also  i  ncorporates  a 
Boltzmann  factor  that  accounts  for  the  thermal  populations  of  the  rovibrational  states  involved  in 
the  t  ransition  w  hich  i  s  t  he  r  eason  f  or  i  t  be  ing  a  f  unction  of  t  emperature.  W  ith  t  hese  c  hanges 
Beer’s  law  becomes: 


1  =  Iae-SIT)‘l--)N’L  ->  A,  =  S(T)g(v - v„)NtL  (2) 

where  Nt  is  the  total  density  of  the  absorbing  species.  The  lineshape  function  is  normalized  such 
that  i  ntegration  ove  r  f requency  o  r  w  avenumber  yields  a  v  alue  of  uni  ty.  T  he  d  ata  h  ere  c  an  b  e 
described  by  a  Voigt  lineshape  function  which  takes  into  account  the  collisional  broadening  and 
the  Doppler  broadening  through  the  convolution  of  a  Gaussian  with  a  Lorentzian.7  It  is  computed 
through  the  use  of  the  Humlicek  algorithm.  This  means  that  if  Equation  2  i  s  integrated  over 
frequency  then  the  frequency  integrated  absorption  is  related  to  the  line  intensity  and  the  column 
density  (NtL)  through  the  following  relation: 


|  Avdv  =  A  =  S(T)NtL  =  aL  (3) 


where  t  he  frequency  integrated  ab  sorbance,  A,  a  nd  t  he  f  requency  i  ntegrated  a  bsorption 
coefficient,  a  ,  have  been  introduced.  The  frequency  integrated  absorption  coefficient  has  been 
introduced  be  cause  i  t  m  akes  de  scription  of  i  nhomogeneous  flow  fields  p  articularly  s  imple.  In 
that  case,  the  frequency  integrated  absorbance  is  given  by: 

A=  |  a(l)dl  (4) 

path 


This  is  an  innocuously  looking  formula  that  can  actually  be  very  difficult  to  compute  because, 
with  the  adoption  of  the  HITRAN  line  strength  in  Equation  3,  a  ,  is  a  complicated  function  of 
temperature  and  pressure.  F  or  example,  if  significant  variations  in  the  lineshape  occur  because 
of  variations  in  absorbing  species  temperature,  density,  and  total  density  then  computation  of  A 
will  require  more  than  a  single  Voigt  lineshape  to  fit  the  absorption  peak.13,14  This  is  found  not  to 
be  problem  when  examining  the  isolator  data  based  on  the  fit  residuals.  Another  obvious  aspect 
is  th  at  if  lo  cal  flow  in  formation  is  d  esired  th  en  it  b  ecomes  n  ecessary  to  u  se  s  ome  f  onn  o  f 
tomographic  i  nversion  or  f  low  f  ield  m  odel  t  o  he  lp  detennine  a.13’15'19  For  t  he  i  solator 
experiments,  CFD  simulations  are  used  to  compute  the  frequency  integrated  absorption  that  can 
then  be  used  to  compare  the  simulations  with  the  measurements.  This  is  done  in  two  steps:  (1) 
an  ima  ge  o  r  m  ap  o  f  t  he  f  requency  i  ntegrated  absorption  c  oefficient  i  s  ge  nerated  f  rom  t  he 
temperature  and  densities  computed  by  CFD,  and  (2)  the  path  i ntegration  i s  done  numerically 
using  a  t  rapezoidal  i  ntegration  r  outine.  T  hen  the  f  requency  i  ntegrated  a  bsorption  for  e  ach 
absorption  feature  can  then  be  input  into  the  routine  for  computing  the  temperature  and  species 
concentration  from  a  Boltzmann  plot.  This  allows  direct  comparisons  between  the  CFD  and  path 
integrated  measurements  to  be  made. 
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3.2  Spectral  Line  Parameters  and  Sample  Spectra 


Water  va  por  ha  s  s  everal  s  trong  r  ovibrational  s  pectroscopic  t  ransitions  f  rom  t  he  vi  sible  t  o  t  he 
mid-IR.  The  availability  of  compact,  spectroscopic  quality  lasers  at  telecom  wavelengths,  1250- 
1650  nrn  ,  ha  ve  m  ade  t  he  2  vi,  2  V3,  a  nd  V1+V3  vibrational  b  ands  p  articularly  attractive  f  or 
diagnostic  applications.  The  abundance  of  water  vapor  as  a  primary  combustion  product  and  the 
relatively  large  line  strengths  for  these  transitions  enable  these  transitions  to  be  used  with  good 
sensitivity  for  many  applications. 

Single  mode,  fiber-coupled,  DFB  lasers  are  used  to  generate  near-infrared  radiation  to  measure 
H20  line  intensities  in  the  spectral  regions  7184.2-7186.4  cm'1,  7159.5-7161.5  cm'1,  and  7179.5- 
7181.8  cm'1.  Principally,  P  and  Q  branch  transitions  of  the  vi  +V3  band  have  been  used  in  this 
study  as  identified  by  the  lowercase  letters  in  figure  1 .  The  data  was  taken  during  time  when  the 
isolator  was  not  back  pressured.  As  can  be  seen  in  the  figure,  peaks  a-d  are  quite  strong  and  they 
are  used  to  detennine  the  temperature  and  water  concentration  in  the  isolator.  Peak  a  is  also  used 
to  determine  the  static  pressure  in  the  isolator.  An  additional  side  peak  is  close  to  feature  c  but 
the  one  at  7181.15  cm'1  is  much  stronger.  Additionally,  this  side  peak  shares  the  same  initial 
state  with  peak  d  so  its  area  can  easily  be  subtracted  from  the  peak  integration.  Finally,  peaks  e-f 
have  much  higher  initial  state  energies  and  so  are  only  used  for  much  higher  temperatures.  The 
data  quality  of  the  time-averaged  data  is  excellent  as  shown  in  Figure  1  with  the  residuals  being 
less  than  1%  of  the  maximum  peak  amplitude  when  using  a  multi-peak  fit  Voigt  algorithm.  Here 
the  r  aw  d  ata  is  in  itially  d  ivided  b  y  th  e  r  eference  beam  m  easurement  after  b  oth  h  ave  b  een 
appropriately  s  caled  1 0  match  i  n  t  he  regions  ou  tside  of  t  he  absorption  peaks.  T  hen  a  c  ubic 
polynomial  plus  Voigt  peaks  are  fit  to  the  resulting  spectrum.  The  cubic  polynomial  corrects  for 
any  nonlinearity  that  is  present  as  the  laser  rapidly  scans  across  the  spectral  region. 
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Figure  1:  Time  averaged  (200  averages)  spectra,  fit,  and  residuals  for  the  three  spectral 

regions  during  run  AC  (see  Table  3) 
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The  line  parameters  are  taken  from  the  HITRAN  2004  database  and  are  listed  in  T able  1 .  The 
line  i  ntensities  a  t  d  ifferent  t  emperatures  a  re  de  rived  us  ing  t  he  s  ame  pr  ocedures  a  nd  pa  rtition 
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function  expressions  as  defined  in  HITRAN.1 1  Here  the  line  parameters  have  been  taken  from 
the  H ITRAN  2004  da  tabase  be  cause  i  t  i  s  t  he  most  c  omplete  s  ource  a  nd  i  t  c  ompares  w  ell  t  o 
recent  s  tudies.  ’  Some  of  t  hese  t  ransitions  ha  ve  be  en  s  tudied  b  y  s  everal  a  uthors,  a  nd  t  he 
variation  among  them  is  between  2-5%  depending  on  t  he  parameter.  ’  '  In  this  study,  the 
HITRAN  parameters  are  used  which  means  systematic  errors  less  than  5%  may  be  present  in  the 
following  analysis  based  on  these  line  parameters.21  It  is  possible  to  reduce  the  systematic  error 
through  careful  self-calibration  of  the  system  over  the  temperature  and  pressure  range  of  interest, 
however,  this  was  not  done  here  because  it  requires  a  special  setup. 


Table  1:  Spectroscopic  Parameters  of  Water  Lines  Used  in  Study 


k  (nm) 

Vio  (cm1) 

S  at  296  K 

(cm'Vmolecule- 

cm'2) 

E” 

(cm1) 

[Vi’,v2’,v3’] 

J’Ka’Kc’ 

J”Ka”Kc” 

Label 

1391.67 

7185.60 

7.947e-22 

1045.06 

1  0  1 

660 

66  1 

a 

1396.37 

7161.41 

1.174e-20 

224.84 

1  0  1 

3  1  3 

4  1  4 

b 

1392.53 

7181.16 

1.505e-20 

136.76 

1  0  1 

2  02 

3  03 

c 

1392.81 

7179.75 

2.299e-22 

1216.19 

1  0  1 

762 

76  1 

d 

1396.49 

7160.81 

5.406e-24 

1874.97 

1  0  1 

1064 

1065 

e 

1392.64 

7180.61 

3.038e-23 

1477.30 

1  0  1 

955 

954 

f 

1392.55 

7181.10 

6.858e-23 

1216.19 

1  20 

770 

76  1 

c’ 

3.3  Determination  of  Temperature  from  a  Boltzmann  Plot 


The  determination  of  path  averaged  static  temperature  is  accomplished  by  measuring  multiple 
absorption  features  of  the  same  molecule.  The  frequency  integrated  absorption  is  normalized  by 
the  cross  section  and  degeneracy  factor  for  each  transition  and  then  plotted  versus  the  initial  state 
energy  di  vided  b  y  k  B.  In  te  nns  o  f  th  e  lin  e  in  tensity  a  bove,  th  is  r  esults  in  th  e  f  ollowing 
relationship: 
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(5) 


where  Q(T)  is  the  partition  function  for  the  absorbing  species  and  E”  is  the  initial  state  energy. 
At  the  temperatures  and  transition  frequencies  studied  here  the  stimulated  emission  tenns  can  be 
ignored  (i.e.  e~hv">IKT  «  0 ).  F  igure  2  shows  a  typical  Boltzmann  plot  for  the  Mach  2  flow  (no 
shock)  taken  during  run  condition  AB.  Here  four  water  transitions  are  used  to  compute  the  path 
integrated  t  emperature.  T  he  cu  rve  s  hows  ex  cellent  1  inearity  which  i  s  i  ndicative  o  f  a  n  early 
uniform  temperature  distribution.  The  slope  of  the  curve  is  inversely  proportional  to  temperature 
and  can  be  directly  determined  from  the  linear  fit.  In  contrast,  the  partition  function1 1  must  be 
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used  to  determine  the  species  concentration  either  from  the  slope  of  the  curve  or  through  using 
the  measured  temperature  and  the  frequency  integrated  absorption  of  one  of  the  transitions. 


Figure  2:  Boltzmann  plot  for  Run  AC  using  the  LOS  average  data  from  a  single  beam 
path.  The  error  bars  are  standard  errors  derived  from  the  curve  fitting  procedure  and  do 
not  include  errors  associated  with  spatial  variations  in  the  flow  as  well  as  the  effect  of 

systematic  errors  in  the  line  parameters. 


3.4  Determination  of  Pressure  from  Line  Width 


Static  pressure  can  be  determined  from  the  collisional  broadening  of  the  absorption  line.  The 
width  of  the  peak  is  proportional  to  two  factors.  T  he  first  factor  is  Doppler  broadening  which 
occurs  b  ecause  t  he  t  hermal  di  stribution  of  t  he  ve  locity  o  f  t  he  m  olecules  i  n  t  he  gas.  This  is 
independent  of  the  states  involved  in  the  transition  and  takes  the  form  of  a  Gaussian  lineshape 
with  a  half  width  at  half  maximum  (HWHM)  of: 


Yd  =3.5812xlO~7v0  (6) 

V  M 

where  Vo,  the  transistion  frequency,  is  given  in  wavenumbers,  T  is  the  gas  temperature  in  K,  and 
M  is  the  molecular  mass  in  grams/mole.  A  t  the  temperatures  in  the  isolator  (400-650  K),  yo 
ranges  from  0.012  -  0.015  cm'1  for  lines  such  as  the  7185.6  cm'1  feature  which  is  of  the  same 
magnitude  as  the  collisional  broadening.  It  is  therefore  desirable  to  first  estimate  the  temperature 
using  the  Boltzmann  plot  and  fix  yD.  This  allows  the  determination  of  the  collisional  broadening 
from  th  e  Voigt  fit  w  ith  a  fixed  D  oppler  br  oadening  w  idth  de  tennined  f  orm  t  he  B  oltzmann 
analysis. 
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Figure  3.  (a)  2D  CFD  temperature  profile  in  the  isolator  (flow  direction  is  out  of  the 
page)  for  run  condition  AC  as  well  as  the  LOS  that  is  being  used  to  calculate  the  path 
averaged  measurement  (b)  CFD  simulated  lineshape,  curvefit  to  a  single  Voigt  profile 
and  residuals,  and  the  LOS  pressure  and  temperature  from  the  fit  compared  to  the 
pressure  and  temperature  at  the  centerline  of  the  isolator. 


An  i  mportant  c  onsideration  w  hen  de  termining  p  ressure  i  n  t  he  m  anner  d  escribed  above  i  s  t  he 
approximation  of  the  absorption  line  by  a  single  V oigt  peak.  In  reality,  in  the  case  of  Mach  2 
duct  flow  (no  shock)  the  temperature  and  density  change  near  the  sidewall  due  to  the  boundary 
layer  with  only  slight  variations  occurring  in  pressure  (on  the  order  of  1%).  In  the  presence  of 
shock  structure,  it  is  necessary  to  validate  the  appropriateness  of  using  a  single  Voigt  profile  with 
significant  p ressure,  temperature,  and  density  inhomogenities  along  the  LOS.  In  severe  cases, 
the  lineshape  should  be  represented  as  a  superposition  of  Voigt  profiles.  To  assess  the  effect  of 
the  approximation  here,  simulations  have  been  conducted  using  the  results  of  CFD  (Reference  6 
addresses  how  the  CFD  flow  simulations  were  calculated).  In  figure  3,  the  upper  plot  shows  a 
pseudo-color  plot  of  the  temperature  for  a  vertical  cross  section  of  the  isolator  duct.  As  can  be 
seen  here,  there  is  a  t  emperature  gradient  of  approximately  200  K  through  the  boundary  layer. 
Because  of  this  temperature  gradient,  the  path  integrated  temperature  is  higher  (630  K)  than  the 
temperature  at  the  centerline  of  the  duct  which  is  6 13  K.  Despite  the  temperature  gradient  and 
the  subsequent  density  gradient  (pressure  gradients  are  small  here  ~  1%),  the  single  line  V oigt 
does  an  excellent  job  fitting  the  composite  lineshape  as  indicated  by  the  residuals  in  this  figure. 
The  oscillations  in  the  residuals  is  the  result  of  trying  to  fit  a  composite  lineshape  with  a  single 
Voigt  as  described  elsewhere.14  Perhaps,  the  most  s  ignificant  aspect  of  the  plot  though  is  the 
excellent  agreement  between  the  static  pressure  determined  from  the  fit  and  the  static  pressure  at 
the  d  uct  cen  terline.  This  s  hows  t  hat  the  bounda  ry  1  ayer  e  ffects  on  t  he  c  urrent  s  tatic  pr  essure 
measurements  should  be  minimal  for  the  Mach  2  flow. 
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Collisional  broadening  is  dependent  on  t  he  gas  molecules  involved  in  the  collisions  as  well  as 
their  quantum  states.  This  dependence  necessitates  the  use  of  different  Lorentzian  broadening 
parameters  f  or  each  t  ransition,  a  nd  t  hese  pa  rameters  m  ust  be  de  tennined  f  or  e  ach  t  ype  of 
collision  partner  over  the  relevant  range  of  conditions  for  the  gas  mixture  being  probed.  A  molar 
based  addition  of  the  broadening  coefficients  for  a  given  temperature  yields  the  total  collisional 
broadening  co  efficient,  yL.  F  or  th  e  v  itiated  gas  mixture  p  resent  in  th  e  i  solator,  th  e  c  ollisional 
broadening  coefficient,  yL,  is  defined  in  equation  7  as 
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(7) 


where  x  is  t  he  m  ole  f  raction  o  f  t  he  p  articular  s  pecies  an  d  t  he  s  pecies  s  pecific  b  roadening 
coefficients  for  the  major  flow  constituents  are  used,  Tref  is  the  reference  temperature  (296  K 
here),  T  is  the  gas  temperature,  and  n  is  the  temperature  coefficient.  The  particular  mole  fraction 
of  th  e  ma  jor  f  low  c  onstituents  is  d  etermined  f  rom  c  alibrated  ma  ss  f  low  a  nd  th  e  me  asured 
>99.9%  c  ombustion  e  fficiency  o  f  t  he  burner.  B  ecause  f  our  br  oadening  c  oefficients  a  nd 
temperature  coefficients  are  required,  it  is  difficult  to  find  a  single  literature  source  that  has  a 
complete  set.  The  computational  work  of  Delaye  et.  al24  lists  the  required  quantities.  However, 
as  pointed  out  by  Nagali  et.  al'  the  temperature  coefficients  can  have  significant  error  for  the 
7185.6  cm1  line.  T  he  self-broadening  c oefficient  of  H 2O  found  in  the  w ork  of  N agali  e t  al 
disagrees  with  a  m  ore  recent  measurement  of  Lepere,  et  al.22  In  addition,  neither  experimental 
author  r  eports  a  t  emperature  co  efficient  f  or  t  he  0  xygen  1  ine  b  roadening  p  arameter.  T  hese 
disagreements  t  hough  a  re  not  unc  oinmon  a  nd  a  c  ritical  a  lthough  s  omewhat  da  ted  r  eview  of 
problems  associated  with  collisional  broadening  parameters  for  water  vapor  lines  is  given  in  the 
work  of  Gamache,  et  al.26  As  stated  therein  disagreements  within  line  parameters  typically  can 

lead  t  0  va  nations  i  n  e  stimated  pr  essure  on  t  he  or  der  of  10%  .  H  ere  the  m  ost  c  urrent  1  ine 
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broadening  parameters  from  the  literature  sources  cited  above  and  a  recent  study  by  Durry,  et  al 
of  N2  and  O2  broadening  are  combined.  The  line  parameters  for  the  7185.6  cm'1  line  used  in  this 
study  are  listed  in  Table  2. 


Table  2:  Broadening  Parameters  and  Temperature  Coefficients  for  line  a  (7185.6  cm'1) 


Perturbing  Species 

y  at  296  K  (cm'1  atm'1) 

Temperature  Coefficient  n 

H20 

0.23  31 

0.56  33 

n2 

0.0456  18 

0.64  33 

02 

0.0271  18 

0.49  17 

co2 

0.086  33 

0.74  33 
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4.0  DISCUSSION  OF  RESULTS 


4.1  Measurements  in  an  Isolator 

In  this  study,  the  static  pressure,  temperature,  and  density  fluctuations  in  a  supersonic  shock  train 
of  a  model  scramjet  isolator  is  explored  using  a  TDLAS  sensor  based  on  16  LOS.  This  particular 
problem  has  been  chosen  for  several  reasons.  Significant  interest  exists  in  the  development  of  air 
mass-flux  orm  ass-flow  s  ensors  fo  r  i  n-flight  m  easurements  ba  sed  on  T  DLAS  t  echnology  f  or 
aerospace  applications.  Although  excellent  measurement  accuracy  to  the  level  of  1%  has  been 
previously  r  eported  f  or  1  ow  s  peed  f  lows,  t  he  c  onditions  f  ound  i  n  s  cramjet  i  solators  pos  e 
significant  ne  w  c  hallenges.  T  he  f  low  i  s  t  ypically  t  urbulent  be  cause  of  t  he  hi  gh  R  eynolds 
numbers  a  nd  s  hock  s  tructures  can  b  e  p  resent  i  n  t  he  flow  ba  sed  on  t  he  e  ngine  ope  rating 
conditions  in  addition  to  the  problems  of  boundary  layer  effects.  By  combining  14  LOS  with  2 
diagonal  LOS  for  velocity  measurements  it  is  possible  to  begin  to  develop  an  understanding  of 
the  influence  shock  trains  on  mass-flux  measurements  although  this  will  not  be  the  current  focus 
of  the  two  part  series.  Instead  the  interest  is  on  the  structure  of  the  shock  train  itself.  Much  is 
known  about  supersonic  shock  trains  in  ducts  which  result  from  the  interaction  of  shocks  with 
the  boundary  layer.  In  a  supersonic  shock  train,  a  sequence  of  shocks  forms  in  the  streamwise 
direction  t  hat  de  pends  t  o  a  1  arge  e  xtent  on  t  he  M  ach  num  her  a  nd  bou  ndary  1  ayer  t  hickness. 
Although  t  he  f  low  i  s  t  ruly  t  hree  di  mensional  i  n  na  ture,  t  he  dom  inant  va  nations  oc  cur  i  n 
streamwise  d  irection  a  nd  s  o  s  ignificant  in  formation  a  bout  t  he  s  hock  t  rain  s  tructure  can  be 
obtained  t  hrough  pr  obing  t  hat  di  rection  o  f  t  he  flow.  T  his  s  implified  t  he  e  xperimental  s  etup 
considerably  a  nd  a  llowed  t  ests  t  o  be  c  onducted  w  ith  e  xisting  opt  ical  a  ccess  ha  rdware.  In  t  he 
future,  it  should  be  possible  to  extend  these  measurements  by  probing  the  transverse  direction  as 
well  a  1 1  he  c  ost  of  a  m  ore  c  omplex  s  etup.  M  uch  i  nfonnation  ha  s  be  en  obt  ained  us  ing  w  all 
measurements  of  this  phenomena  but  is  limited  to  determining  the  shock  front  location  and  much 
less  information  exists  about  shock  train  inflow  properties  as  well  shock  train  dynamics.  H  ere 
the  goal  is  to  obtain  time  dependent  inflow  data  ( 1  kHz)  at  several  spatial  locations  along  the 
flow  direction  that  can  be  compared  to  theory  as  well  as  computational  fluid  dynamics  (CFD) 
simulations.  T  raditionally  these  flow  properties  have  been  explored  using  particle  velocimetry 
techniques,  high-speed  Schlieren,  and  in-stream  pressure  probes.  These  studies  though  focused 
on  slower  flows  (M-1.6)  than  the  ones  studied  here  (M-2.2)  where  the  shock  train  unsteadiness 
was  much  less  and  the  shock  train  was  known  to  be  a  normal  shock  train.  T  he  use  of  TDLAS 
can  complement  t  hese  measurements,  an  d  i  t  o  ffers  t  he  s  ignificant  ad  vantage  t  hat  i  t  c  an  b  e 
miniaturized  so  that  in-flight  measurements  could  be  possible  in  addition  to  ground  testing. 

In  this  re  port,  t  he  s  ystem  de  sign,  m  easurement  t  echnique  a  nd  M  ach  2  flow  r  esults  w  ill  be 
presented.  This  allows  comparison  of  system  perfonnance  to  previous  measurements  that  have 
been  made  in  the  literature  as  well  as  to  conventional  pressure  transducers  and  thermocouples. 
Also  i  ntroduced  he  re  i  s  t  he  us  e  of  c  ollisional  b  roadening  tom  easure  t  he  s  tatic  pr  essure  i  n  a 
supersonic  flow.  A  lthough  collisional  broadening  has  been  well  known  since  the  beginning  of 
spectroscopic  t  echniques,  t  o  our  know  ledge  i  t  h  as  not  be  en  e  xploited  i  n  t  he  m  easurement  of 
supersonic  flows.  Probably  the  most  significant  reason  for  this  is  that  previous  campaigns  have 
not  been  interested  in  flows  with  substantial  pressure  variations.  However,  there  is  also  the  issue 
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of  a  ccurately  know  ing  t  he  1  ine  br  oadening  p  arameters  w  hich  a  re  continually  i  mproving  w  ith 
more  accurate  measurements  and  calculations  being  made. 1 1  Here  it  will  be  shown  that  accurate 
measurements  of  pressure  can  be  made  if  the  line  parameters  are  well  known. 


4.1.1  Optical  Layout  and  Control/Acquisition  Electronics 


Figure  5  shows  a  schematic  of  the  common  optical  components  used  in  each  experimental  setup. 
Three  di  stributed  f  eedback  di  ode  lasers  (  NEL  NLK1E5GAAA)  al  1  o  perating  i  n  a  d  ifferent 
spectral  region  are  time-multiplexed  in  a  single  period  lasting  1  millisecond  and  combined  using 
a  lx  4  c  ombiner  ( Newport  F  -CPL-B 14350).  T  he  c  ombined  be  am  i  s  t  hen  s  plit  us  ing  a  10/90 
splitter  (Newport  F-CPL-B 12351)  with  the  10%  beam  going  to  a  50/50  splitter  (Newport  F-CPL- 
B 12355)  whose  output  is  routed  to  either  a  SiCfi  (Los  Gatos  Research)  etalon  with  a  2.000+0.002 
GHz  @  1550  nm  free  spectral  range  that  monitors  the  frequency  spectrum  of  the  lasers  or  to  a 
photodetector  (Thorlabs  PDA  400)  that  is  used  as  a  reference.  The  remaining  90%  of  the  beam 
enters  another  50/50  splitter  that  then  connects  to  two  1x8  fiber  optic  splitters  (Newport  F-CPL- 
lx8-OPT)  that  are  then  collimated  (Thorlabs  F240APC-C  8mm  FL).  The  signal  is  attenuated  by 
a  factor  of  -18.5  dB  as  it  passes  through  the  splitters  and  ultimately  each  LOS  will  see  an  average 
power  of  280  pW  if  the  average  input  power  of  the  lasers  is  20  mW.  T  his  power  level  allows 
excellent  signal  to  noise  with  the  detector  arrays  as  can  be  seen  in  figure  1 .  The  collimators  are 
mounted  on  a  pr  ecision  dr  died  ba  seplate  toe  nable  t  ight  be  am  pi  acement.  Because  of  t  he 
physical  limitations  of  the  optical  mounts,  the  spacing  of  the  collimator  arrays  is  slightly  larger 
than  the  3”  photodiode  array  used  for  detection. 

The  raw  transmitted  laser  power  is  measured  simultaneously  on  a  11  16  pa  ths  through  the  use  of 
in-house  custom  designed  InGaAs  photodiode  arrays.  These  custom  designed  arrays  amplify  the 
signal  f  rom  2m  m  di  ameter  FGA21  phot  odiodes  (  FGA  21  f  rom  T  horLabs)  us  ing  a  T  exas 
Instruments  OPA380  transimpedance  amplifier/op  amp  for  each  photodiode.  This  amplifier  with 
90  MHz  of  gain  bandwidth  is  a  single  supply  (+5  V)  and  therefore  requires  a  bias  be  applied  to 
the  photodiode  so  that  input  voltages  are  not  driving  to  0  V  .  The  simplified  schematic  of  the 
optical  amplifier  circuit  is  shown  in  Figure  5.  Several  design  tradeoffs  were  made  in  the  circuit 
design  to  allow  for  an  easily  scalable  system  as  well  as  ease  of  use.  A  256  tap,  200  kQ  digital 
potentiometer  (AD5262)  allows  for  the  gain  to  be  controlled  digitally.  The  leads  to  an  additional 
25  pF  of  parasitic  capacitance  that  decreases  the  bandwidth  and  combined  with  the  photodiode 
junction  c  apacitance  of  500  pF  a  t  0  bi  as  vol  tage  ( at  0.5  V  bi  as  i  t  i  s  1  ess)  yields  a  t  otal  i  nput 
capacitance  of525pF.  T  his  cap  acitance  w  ith  g  ain  s  etting  o  f  3  0  k  Q  results  i  n  an  o  verall 
bandwidth  of  approximately  950  kHz  as  determined  from  the  manufacturer’s  datasheet.  At  200 
kQ  the  bandwidth  of  the  circuit  will  decrease  to  370  kH  z.  A  lthough  this  seems  considerably 
faster  than  the  1  kHz  scan  rate  the  primary  issue  is  the  scan  time  across  the  peak  which  here  is 
approximately  5.6  pS.  The  frequency  response  must  be  faster  in  order  to  get  a  correct  peak  shape 
and  area  which  is  the  case  here  for  the  nominal  30  kQ  feedback  resistor  used;  however,  at  high 
gains  this  can  become  a  problem.  T his  can  be  improved  by  moving  to  faster  amplifiers  with  a 
larger  gain-bandwidth  product.  T  he  use  of  smaller  photodiodes  are  possible  as  well,  however, 
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To  Array  1 


To  Array  2 


Figure  4:  Optical  layout  of  multiplexed  diode  laser  system  and  time  multiplexed  signal 
(inset)  with  peaks  labeled  as  indicated  in  Figure  1  and  Table  1 

noise  due  to  optical  turbulence  increases  as  the  photodiode  area  decreases  if  the  beam  overfills 
the  detector.  S  o,  clearly  improvements  can  be  made  to  these  photodiode  arrays,  however,  the 
signal  is  still  excellent  for  the  size  of  the  measured  absorption  peaks  in  Figure  1 . 


AD5262  0-200  kfi 


Figure  5:  Simplified  schematic  of  the  optical  amplifier  circuit 

The  signal  is  then  digitized  into  14  bits  using  a  NI  PXI-6133  (32  MSamples  onboard  memory)  S 
Series  Multifunction  Digital  Acquisition  (DAQ)  Device  sampling  at  a  rate  of  2.5  MSamples/s. 
The  amplifier  output  as  well  as  all  control  signals  and  power  are  provided  through  a  68  w  ire 
cable  that  connects  directly  to  the  photodiode  array  using  a  SCSI-3  type  connector.  This  along 
with  a  daptable  s  oftware  w  ritten  b  y  u  s  a  Hows  f  or  e  asy  s  calability  o  f  th  e  s  ystem  b  ecause  each 
digitizer  card  and  photodiode  array  provide  an  additional  8  c  hannels.  A  PX1  chassis  is  used  to 
hold  all  of  the  data  acquisition  equipment;  i.e.  computer/controller,  DAQ  cards,  and  the  function 
generators  used  to  drive  the  lasers.  Because  of  the  unsteady  nature  of  the  shock,  a  NI  Intel  Core 
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Duo  2.0  GHz  Processor  with  2GB  of  DDR2  host  RAM  is  used  to  average  200  spectra,  determine 
peak  areas  through  trapezoidal  integration,  and  then  create  a  real  time  (1-5  Hz)  LOS  image  along 
the  flow  direction.  This  provided  feedback  as  to  the  success  of  the  data  runs.  It  also  provides 
control  of  the  current  ramp  of  laser  diodes  through  a  N ational  Instruments  PXI-5402  function 
generator.  T  he  digitizers  and  function  generators  are  synchronized  to  ~  5  nanoseconds  using  a 
10  MHz  backplane  clock  signal.  Because  all  the  channels  are  all  simultaneously  sampled  at  the 
same  time  (18x2.5x2  =  90  MB/s),  the  PCI  bus  in  the  PXI  chassis  cannot  sustain  continuous  data 
transfer  to  the  hard  drive  in  the  current  configuration.  Instead  a  short  segment  of  data  (0.2  s)  is 
captured,  processed,  and  saved  to  the  hard  drive.  This  process  takes  approximately  5  s  and  then 
another  0.2  s  of  data  is  captured,  processed,  and  saved.  A  lthough  the  shock  is  unsteady  it  has 
been  found  that  pushing  the  shock  slowly  from  the  rear  of  the  isolator  to  the  front  enabled  the 
key  shock  features  to  be  measured.  The  data  acquisition/processing  software  was  developed  in 
Lab  Windows  to  allow  a  scalable  architecture  for  developing  TDLAS  based  tomography.  The 
data  a  cquisition  r  ates  w  ill  va  ry  ba  sed  on  t  he  c  hannel  count,  how  ever,  at  1  east  2  s  econds  o  f 
continuous  data  can  be  acquired  because  of  the  large  on-board  memory  of  the  digitizer  cards. 


4.1.2  Isolator  Test  Setup 


Figure  6.  (a)  Schematic  of  the  supersonic  combustion  test  facility  and  (b)  photograph  of  the 
facility  that  indicate  1  -  burner,  2  -  facility  nozzle,  3  -  TDLAS  hardware,  4  -  isolator,  5  - 
extension,  6  -  backpressure  valve,  and  7  -  exhaust  line 

AFRL/RZA  owns  and  operates  two  direct-connect  supersonic  combustion  facilities.  Figure  6 
shows  a  schematic  of  the  Research  Cell  18  where  the  experiments  were  performed.  This  facility 
was  designed  to  allow  basic  studies  of  supersonic  reacting  flows  using  conventional  and  non- 
intrusive  di  agnostic  t  echniques.  It  consists  of  a  natural  gas  fueled  burner,  t  ransition  flange, 
interchangeable  facility  nozzle  (Mach-1. 8  and  2.2  currently  available),  modular  isolator,  modular 
combustor,  instrument  section  and  calorimeter.  The  rig  is  mounted  to  a  thrust  stand 
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capable  of  measuring  thrust  of  up  to  2000  lbf.  The  test  rig  is  provided  with  continuous  flow  of  up 
to  30  lbm/sec  at  750  psia  at  a  temperature  of  up  to  922F,  and  exhausts  into  a  3.5  psia  continuous 
flow  exhauster  system.  The  facility  can  simulate  flight  conditions  from  Mach  3.5  to  Mach  5  at 
flight  dynamic  pressures  of  up  to  2000  p  sf.  F  or  the  current  study,  the  combustor  was  removed 
and  r  eplaced  w  ith  a  1  arge  valve  t  hat  w  as  u  sed  t  o  s  imulate  t  he  p  ressure  r  ise  as  sociated  w  ith 
combustion.  The  isolator  is  fitted  with  1”  high  x  3.4”  wide  optical  access  ports  designed  to  hold 
1.7  cm  (0.68  in)  inch  thick  quartz  windows  or  steel  inserts.  The  quartz  windows  are  wedged  at 
5°  on  the  outer  surface  to  prevent  etaloning.  The  nozzle  and  isolator  are  instrumented  with  static 
pressure  t  aps  and  t  hermocouple  w  ith  a  frequency  r  esponse  of  a  pproximately  10  H  z.  N  o 
difference  in  isolator  perfonnance  has  been  noted  between  the  using  the  quartz  windows  or  the 
thennal  barrier  coated  (TBC)  steel  inserts  despite  the  differences  in  surface  roughness  between 
the  qua  rtz  a  nd  t  he  T  BC .  A  ke  y  e  nabling  a  spect  pr  ovided  b  y  t  his  f  acility  is  t  hat  t  he  vi  tiation 
which  is  used  to  increase  the  flow  enthalpy  to  the  correct  level  provides  the  water  vapor  being 
probed  in  this  experiment.  A  lthough  it  is  possible  to  measure  oxygen,  it  is  difficult  to  achieve 
the  e  xcellent  s  ignal  t  o  noise  r  atio  t  hat  c  an  be  obtained  w  ith  w  ater  a  bsorption  f  eatures.  This 
facility  has  been  described  in  detail  elsewhere.  ’ 


Table  3  Average  Run  Conditions  -  Facility  Nozzle  Mach  Number  =  2.2 


Start  Time 

Incoming  Air 
Temperature  (K) 

Total  Temperature 
(K) 

Total  Pressure 
(psia) 

Water  mole 
percent  (%)* 

Reynolds 

Number 

19:23:07 

300 

667 

50 

3.03 

9.1xl05 

19:36:21 

300 

1016 

50 

6.93 

4.6xl05 

21:26:20 

589 

667 

50 

1.11 

9.1xl05 

21:34:58 

589 

1016 

50 

4.94 

4.6xl05 

*Assumes  dry  facility  air 


In  this  report,  the  analysis  will  focus  on  four  sets  of  data  labeled  runs  AB,  AC,  AF,  and  AG.  The 
burner  settings  as  w  ell  as  t  he  noz  zle  M  ach  nu  mber,  w  ater  m  olar  pe  rcent,  a  nd  t  he  R  eynolds 
number  for  each  case  are  shown  in  Table  3.  The  Reynolds  number  has  been  computed  using  the 
equivalent  duct  diameter  (5.75  cm)  for  a  rectangular  duct  with  cross  section  of  3.81  x  10.16  cm. 
The  amount  of  vitiation  used  for  each  run  depended  on  the  desired  total  temperature  and  whether 
the  air  was  preheated  using  electric  heaters.  T he  temperature  and  total  pressure  showed  only  a 
few  de  grees  a  nd  a  bout  a  ha  if  a  ps  i  va  riation  ove  r  t  he  c  ourse  of  a  da  ta  c  ollection  c  ycle  w  hich 
means  t  hat  t  hey  a  re  a  ccurate  t  o  1%  ove  r  t  he  c  ourse  of  t  he  t  est  r  un.  T  he  w  ater  m  olar 
concentration  though  ha  s  t  o  be  computed  f rom  the  m  ass  flows  of  t  he  c  ompressed  na  tural  gas 
(CNG)  a  nd  ox  ygen.  In  t  his  c  ase,  s  everal  i  ssues  e  xist.  F  irst,  C  NG  i  s  not  pur  e  m  ethane  but 
includes  CCF  and  larger  alkanes  such  as  ethane  and  propane.  The  presence  of  CO 2  will  lower  the 
amount  of  water  in  the  vitiated  flow,  but  the  combustion  of  larger  alkanes  will  increase  it.  The 
CNG  used  in  these  experiments  is  composed  of  95.19%  methane,  2.30%  ethane,  1.28%  COo, 
0.65%  nitrogen,  and  0.58%  higher  alkanes  (propane,  butane,  pentane,  etc.)  based  on  an  average 
of  the  gas  composition  reported  by  the  CNG  provider. 
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Wavenumber  (cm’1) 

Figure  7.  Absorption  spectra  taken  before  runs  AF  and  AB  with  no  vitiation  to  assess  the 

effectiveness  of  the  purge  boxes. 

Second,  in  the  time-averaged  case  shown  here,  a  low  frequency  (-0.05  Hz)  1%  oscillation  was 
observed  in  both  mass  flow  rates  of  the  oxygen  and  CNG  controllers  which  is  averaged  out  here 
because  the  data  run  usually  lasted  on  the  order  of  1  V2  minutes  —  2  V2  minutes.  However,  for  the 
comparisons  made  later  a  0  .2  s  data  segment  was  used  and  the  average  water  mole  percent  was 
computed  for  that  particular  time  from  the  mass  flow  controller  data.  F  inally,  variations  in  the 
incoming  a  ir’s  m  oisture  c  ontent  ha  s  be  en  obs  erved  a  s  c  an  be  s  een  i  n  figure  7  w  hich  s  hows 
measurements  made  to  assess  the  effectiveness  of  the  purge  boxes.  Because  the  purge  boxes  do 
such  a  g  ood  j  ob  a  s  e  videnced  b  y  t  he  a  pparent  n  egative  a  bsorption  pe  ak  a  ssociated  w  ith  da  ta 
taken  be  fore  r  un  A  B,  it  i  s  p  ossible  1 0  u  se  t  he  data  1 0  as  sess  t  he  w  ater  m  ole  p  ercent  0  f  t  he 
incoming  air.  Here  the  negative  absorption  peak  is  actually  due  to  a  small  air  path  (-  1  -2mm)  in 
the  reference  signal  that  is  unmatched  on  t  he  measurement  channels  most  likely  due  to  the  air 
gap  between  the  photodiode  and  the  end  of  the  optical  fiber  (a  Thor-Labs  SM1FCA  adapter  plate 
is  used  to  connect  reference  fiber  to  the  PDA  400).  Because  it  is  of  the  order  of  0.01  -0.03% 
water  molar  percent,  we  ignore  it  in  our  calculations.  H  owever,  as  can  be  seen  in  figure  7  t  he 
data  collected  before  run  AF  had  much  greater  water  content.  H  ere  the  water  concentration  of 
2.0  ±  0.1  x  10  molecules/cm  ,  the  static  temperature  of  275  ±  15  K,  and  the  static  pressure  of  4 
±  0.2  psia  indicates  that  the  incoming  air  has  a  water  mole  percent  of  0.28%  ±  0.03%.  Although 
this  is  a  very  small  water  mo  lar  p  ercent,  it  is  very  important  when  assessing  the  ability  of  the 
TDLAS  to  measure  the  water  concentration  in  the  case  of  run  AF  where  this  effect  would  lead  to 
a  25%  error  if  not  corrected. 
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Figure  8.  A  top  view  of  the  experimental  setup  with  an  overlay  depicting  the  shock  train 
structure  in  the  isolator  is  shown.  The  measurement  plane  occurs  at  mid-height  in  the  isolator. 
For  clarity,  only  8  beams  (red  lines)  going  from  the  bottom  purge  box  to  the  top  are  shown. 
However,  another  8  beams  from  the  top  to  the  bottom  box  are  also  used.  The  indicated  flow 
direction  is  what  will  be  used  in  subsequent  figures. 


Figure  8  shows  t  he  ex  perimental  1  ayout  a  s  i  mplemented  dur  ing  t  he  t  est  c  ampaign.  T  he  m  ost 
significant  changes  in  the  flow  structure  occur  in  the  streamwise  direction  and  consist  of  high 
and  low  regions  of  density,  pressure  and  temperature  in  the  shock  train.  It  is  the  goal  here  to 
detect  these  variations  in  good  detail  in  the  streamwise  direction  with  an  approximate  spacing  of 
the  14  LOS  of  5  mm.  Obviously,  higher  aspect  that  can  be  seen  from  the  image  is  the  reason  for 
the  i  ndicated  pr  obing  geometry  of  14  be  ams  be  ing  aligned  ne  arly  p  erpendicular  t  o  t  he  flow 
direction.  The  most  significant  changes  in  the  flow  structure  occur  in  the  streamwise  direction 
and  consist  of  high  and  low  regions  of  density,  pressure  and  temperature  in  the  shock  train.  It  is 
the  goal  h  ere  tod  etect  t  hese  v  ariations  in  good  d  etail  in  th  e  s  treamwise  d  irection  w  ith  a  n 
approximate  s pacing  of  the  14  L  OS  of  5  mm.  Obviously,  hi  gher  r esolution  w ill  i mprove  t he 
quality  of  the  data.  H  owever,  significant  aspects  of  the  shock  structure  can  be  measured.  T  he 
photodiode  arrays  as  w  ell  as  co  llimator  ar  rays  ar  e  m  ounted  o  n  s  tandard  1  2”xl8”  o  ptical 
breadboards  that  are  inserted  into  aluminum  purge  boxes.  These  purge  boxes  were  sealed  except 
for  inlets  and  exit  ports  for  N2  purge  gas  and  a  face  plate  that  had  a  r  ectangular  hole  cut  into 
slightly  larger  than  the  isolator  windows.  The  front  of  each  box  was  press  fit  against  a  ceramic 
firebrick  cutout  that  insulated  the  purge  box  from  the  isolator  wall  that  could  reach  temperatures 
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of  up  t  o  700  -800  K  d  uring  hi  gh-enthalpy  r  uns.  Finally,  a  ma  sk  in  sert  th  at  r  estricted  th  e 
photodiodes  field  of  view  was  found  necessary  for  runs  AC  and  AG  due  to  black-body  radiation 
associated  w  ith  t  he  i  nterior  w  alls  of  t  he  i  solator.  T  he  box  es  w  ere  m  ounted  on  a  m  uch  1  arger 
breadboard  that  was  attached  to  the  thrust  stand  rails  that  allowed  the  whole  setup  to  move  with 
the  test  rig  as  it  thermally  expanded.  Compressed  nitrogen  was  used  to  pressurize  a  LN2  dewar 
whose  blow  off  was  used  to  cool  and  purge  the  boxes  of  ambient  air.  As  mentioned  previously, 
excellent  purge  was  achieved  and  no  evidence  of  water  vapor  due  to  the  regions  outside  the  flow 
was  observed  during  the  experiments. 
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Figure  9  (a)  Boltzmann  plot  for  one  scan  (0.3  ms),  (b)  temporal  variation  of  temperature  versus 
time,  (c)  spectrum  for  one  scan,  and  (d)  pressure  versus  time  all  for  the  LOS  at  x  =  3.43  cm  from 

the  upstream  edge  of  window 
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4.2  Discussion  of  Isolator  Measurements 


The  assessment  of  the  measurement  accuracy  and  precision  are  necessarily  more  difficult  for  the 
case  o  f  t  he  M  ach  2  ( no-shock)  i  solator  flow  t  han  f or  s  teady  flow  pr  oblems  s  uch  a  s  1  aminar 
flames  which  have  been  studied  extensively.  Here  issues  such  as  having  high-quality  data  that 
can  be  used  to  compare  the  measurements  with,  flow  unsteadiness,  and  small  spatial  variations  in 
the  flow  field  all  will  impact  the  accuracy  as  well  as  precision  of  the  measurements.  In  addition, 
systematic  e  rrors  due  toe  rrors  i  n  t  he  s  pectroscopic  pa  rameters  w  ill  a  ffect  t  he  co  mparison. 
Nevertheless,  as  will  be  shown  here  excellent  comparisons  exist  between  static  temperatures  and 
calculations,  s  tatic  p  ressures  an  d  av  erage  w  all  pressure  m  easurements,  an  d  t  he  w  ater  m  ole 
percent  and  complete  combustion  assumptions. 
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Figure  10  Streamwise  variation  of  time  and  path  averaged  static  temperature,  water 
concentration,  and  static  pressure  for  run  conditions  (a)  AB  (b)  AC  (c)  AF,  and  (e)  AG.  Scales  are 
taken  to  be  ±10%  of  the  mean.  The  displayed  error  bars  are  fit  error  bars  and  do  not  reflect 

inaccuracies  due  to  line  parameters. 


Figure  9  shows  the  measurements  for  a  single  LOS  located  3.43  cm  from  the  upstream  edge  of 
the  w indow  dur  ing  r un  AB  ( no  ba ckpressure  c  ondition).  At  a ny  given  time ,  th e  d ata  is  of 
excellent  q  uality  and  t  he  t  emperature  an  d  water  co  ncentration  can  bed  educed  from  t  he 
Boltzmann  pi  ot  a  s  w  ell  a  s  pr  essure  f  rom  t  he  1  inewidth  of  1  ine  a.  T  he  fluctuations  in  s  tatic 
temperature  an  d  p  ressure  are  o  ft  he  o  rder  3%  a  nd  r  effect  v  ariations  in  th  e  rig  o  perating 
conditions  as  well  as  the  baseline  estimation  procedure  used  to  process  the  data. 
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As  can  be  seen  in  Figure  10,  e  xcellent  precision  is  obtained  between  identical  flow  conditions 
despite  significant  time  difference  between  runs  and  the  use  of  unheated  or  preheated  air  to  the 
burner.  Figure  10  also  shows  how  the  time  averaged,  spatially  resolved  data  varies  as  a  function 
of  streamwise  position  in  the  duct.  T  he  maximum  value  of  the  amplitude  variations  is  always 
less  than  10%  of  the  mean  value  (scale  is  ±10%)  and  often  much  less.  The  data  for  this  analysis 
was  acquired  by  first  averaging  the  spectra  from  the  200  di  fferent  scans  in  one  data  segment. 
Multi-peak  Voigt  fits  were  then  used  to  fit  each  spectrum  in  the  manner  described  in  section  II 
and  obt  ain  t  he  da  ta  i  n  f  igure  10.  T  he  da  ta  s  egment  w  as  chosen  bye  xamining  each  r  un  time 
during  the  course  of  each  run  when  it  was  expected  that  no  s  hock  was  present  in  the  window. 
However,  in  case  of  Run  AG,  the  very  front  edge  of  a  shock  is  beginning  to  enter  the  rear  of  the 
window  as  confirmed  by  an  inspection  of  the  isolator  pressure  taps.  All  cases  show  that  the  LOS 
closest  to  the  rear  of  the  window  exhibits  a  rise  in  pressure  compared  to  the  other  LOS.  This  also 
appears  to  be  the  case  for  the  pressure  tap  data  although  it  has  a  much  sparser  resolution.  It  may 
be  the  case  that  a  weak  shock  is  attaching  near  the  rear  of  the  window  where  there  is  a  small 
break  between  the  window  and  the  TBC  coated  sidewall  of  the  isolator.  The  error  bars  in  figure 
10  correspond  to  the  standard  statistical  error  bars  from  a  curve  fit.  These  error  bars  agree  with 
the  results  reported  in  Figure  9  because  the  error  in  the  mean  in  that  case  roughly  corresponds  to 
the  s  tandard  e  rror  ba  rs  i  n  F  igure  10.  N  otice  t  hat  a  lthough  s  tandard  e  rrors  a  re  b  eing  reported 
nothing  i  s  be  ing  s  aid  a  bout  c  onfidence  i  ntervals.  T  his  r  equires  knowledge  of  t  he  unde  rlying 
distribution  of  the  data  which  is  difficult  to  detennine  and  most  likely  not  Gaussian.  However, 
the  s  mall  s  tandard  e  rrors  s  uggest  th  at  s  omething  o  ther  th  an  s  tatistical  e  rrors  are  causing  th  e 
spatial  oscillations  such  as  local  flow  fluctuations.  B  y  averaging  over  these  spatial  oscillations 
an  increase  in  measurement  precision  is  achieved. 


Table  4  Static  Pressure  Comparison 


Run 

TDLAS  (psia) 

Pressure  Taps  (psia) 

%  Difference  [TDLAS-P.  Taps] 

AB 

6.39  ±0.32 

6.59  ±0.20 

-3.0  % 

AC 

6.35  ±0.32 

6.49  ±0.09 

-2.2  % 

AF 

6.55  ±0.33 

6.65  ±0.08 

-1.5  % 

AG 

6.56  ±0.33 

6.58  ±0.14 

-0.3  % 

Figure  1  1  shows  that  spatial  oscillations  also  exist  in  the  sidewall  pressure  measurements.  A  s 
can  be  seen  here  during  run  AC,  the  pressure  fluctuates  along  the  duct.  The  error  bars  indicate 
the  calibration  error  of  the  pressure  transducers  so  the  oscillations  that  are  being  observed  are 
due  t  o  flow  f ield  va  nations.  T  his  s  hould  be  e  xpected  f or  a  t  urbulent  flow.  In  an  e  ffort  t  o 
eliminate  these  fluctuations,  the  three  pressure  taps  that  span  the  window  region  ( 1  bottom  side 
and  2  t  op  side)  are  averaged  together  and  compared  to  the  spatially  and  temporally  averaged 
TDLAS  m  easurements.  T  he  co  mparison  b  etween  t  he  T  DLAS  m  easurements,  p  ressure  t  ap 
measurements,  and  CFD  is  presented  in  Table  4.  As  can  be  seen  in  the  table,  the  agreement  is 
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excellent  with  the  worst  disagreement  being  3%  for  run  AB  and  the  average  disagreement  being 
1.75%  between  TDLAS  and  the  average  pressure  tap  pressure.  This  is  less  than  the  error  bars  for 
the  T  DLAS  m  easurements  w  hich  ha  ve  be  en  t  aken  t  o  be  ±  5%  (  0.32  psia)  due  t  o  pos  sible 
systematic  errors  in  the  line  broadening  parameters  in  accordance  with  reference  34.  The  errors 
for  the  pressure  taps  are  the  standard  deviation  of  the  mean  of  the  three  pressure  taps  located 
within  the  window  region.  Note  that  the  TDLAS  data  shows  remarkable  similarity  between  runs 
AB  &  AC  and  runs  AF  &  AG.  T  he  statistical  precision  of  the  TDLAS  data  can  be  determined 
from  figure  10  to  be  approximately  0.04-0.07  psia  so  it  would  be  much  better  than  the  pressure 
transducer  data  which  only  exhibits  a  precision  of  0.1  -  0.2  psia.  At  the  current  level  of  analysis, 
though,  it  will  be  difficult  to  determine  the  sources  of  these  differences  because  it  could  be  due  to 
a  large  number  of  factors  involved  in  accurately  determining  the  facility  operating  conditions. 

The  s  tatic  t  emperature  o  f  t  he  f  low  can  o  nly  b  e  co  inpared  at  t  he  cu  rrent  t  ime  t  o  t  he  C  FD 
simulations.  In  order  to  do  t  his  properly,  the  path  integrated  temperature  must  be  determined 
from  the  CFD  as  described  in  section  II  part  C. 
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Figure  11.  Pressure  measurements  using  conventional  static  pressure  taps  on  the  bottom  side  of 
the  isolator  as  well  as  the  window  region  near  the  time  the  Mach  2  flow  TDLAS  data  was  taken 

for  run  AC. 


The  c  omparison  i  s  s  hown  i  n  T  able  5  .  A  s  can  b  e  s  een  f  rom  t  his  t  able,  e  xcellent  ag  reement  i  s 
found  between  the  computed  temperatures  and  the  measured  ones  with  a  maximum  difference  of 
3.3%  and  an  average  difference  of  2.3%.  T  he  error  bars  have  been  computed  based  on  a  ±3% 
systematic  error  plus  the  statistical  error  evident  in  figure  10. 

Table  5  a  iso  s  hows  t  he  comparison  be  tween  t  he  f  acility  ba  sed  c  omputations  a  nd  t  he  T  DLAS 
computations  of  water  mole  percent  in  the  vitiated  flow.  The  TDLAS  results  for  the  water  mole 
percent  are  based  off  of  the  TDLAS  static  pressures  in  Table  4,  TDLAS  static  temperatures  in 
Table  5,  a  nd  T  DLAS  w  ater  c  oncentrations.  T  he  i  deal  ga  s  1  aw  i  s  us  ed  t  o  c  ompute  t  he  t  otal 
densities  f  rom  t  he  T  DLAS  s  tatic  p  ressure  and  temperature,  an  d  the  m  easurement  er  rors  ar  e 
propagated  through  using  standard  methods.  T  he  facility  computations  are  based  on  flow  rates 
through  f acility  pr  ovided  m  ass  flow  c  ontrollers  t  hat  pr  ovide  t  he  flow  r  ate  of  a  ir,  c  ompressed 
natural  gas,  and  oxygen  to  the  burner.  The  molar  fraction  of  water  is  then  computed  based  on 
complete  combustion  of  the  CNG  mixture.  Because  the  CNG  is  not  pure  as  detailed  earlier,  the 
maximum  pos  sible  e  rror  of  8%  w  as  be  en  us  ed  t  o  de  tennined  t  he  e  rrors  f  or  t  he  f  acility 
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computations.  I  n  a  ddition,  t  he  m  easured  0.28  %  w  ater  m  ole  pe  rcent  has  be  en  added  i  n  t  he 
computations  f  or  r  uns  AF  a  nd  AG  t  o  a  ccount  f  or  t  he  w  ater  pr  esent  i  n  t  he  f  acility  air.  T  he 
agreement  b  etween  t  he  TDLAS  an  d  t  he  c  alculated  w  ater  m  ole  p  ercent  methods  i  s  q  uite  good 
despite  the  approximations  being  made  with  the  facility  agreeing  with  the  TDLAS  computations 
to  the  level  of  7%  except  for  the  case  of  run  AF.  There  the  difference  in  these  measurements  is 
on  t  he  or  der  of  t  he  f  acility  air  c  orrection  t  o  t  he  s  mall  a  mount  of  w  ater  pr  oduced  unde  r  t  hose 
conditions.  In  general,  the  agreement  is  within  the  error  of  the  facility  computations. 

Before  c  oncluding  t  his  s  ection,  i  t  i  s  us  eful  to  c  ompare  t  he  m  easurements  t  o  pr  evious 
measurement  campaigns  in  similar  facilities.  In  p articular,  a  straightforward  c omparison  w ith 
the  work  of  Upschulte,  et  al.  is  useful  to  see  how  the  measurement  accuracy  and  precision  have 
increased.  F  or  that  particular  case,  two  LOS  were  used  to  probe  the  flow  along  an  18  cm  path 
which  is  almost  a  factor  of  two  longer  than  the  paths  used  in  the  current  experiments  (10.14  cm). 
In  t  hat  w  ork,  the  m  easured  s  tatic  t  emperatures  w  ere  f  ound  t  o  a  gree  w  ith  t  he  p  redicted  s  tatic 
temperatures  w  ith  in  7  -1 1%  w  hich  i  s  a  f  actor  of  2  -3  mo  re  th  an  th  e  ma  ximum  d  isagreements 
observed  he  re.  In  a  ddition,  t  he  w  ater  m  ole  pe  rcent  w  as  obs  erved  t  o  di  ffer  f  rom  t  he  pr  edicted 


Table  5  Path  Averaged  Static  Temperature  Comparison  and  Water  Mole  Percent 


Temperature 

Water  mole  percent 

Run 

TDLAS  (K) 

CFD  (K) 

%  Difference 
[TDLAS-CFD] 

TDLAS  (%) 

Facility 

Computation  (%) 

%  Difference 
[TDLAS-Facility] 

AB 

398  ±4 

389 

2.2  % 

3.30  ±0.20 

3.07  ±0.25 

7.32 

AC 

652  ±  10 

630 

3.3  % 

6.83  ±0.38 

7.09  ±0.57 

-3.7 

AF 

399  ±4 

389 

2.6  % 

1.68  ±0.10 

1.44  ±0.28 
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AG 

638  ±  11 

630 

1.25  % 

5.04  ±0.29 

5.27  ±0.50 

-4.4 

water  f  ractions  by  a  fa  ctor  o  f  1  7-56%  w  hich  i  s  a  f  actor  o  f  1  -3  mo  re  th  an  th  e  ma  ximum 
disagreement  observed  in  this  work.  However,  as  observed  here,  the  lowest  disagreements  were 
observed  with  t  he  hi  ghest  w  ater  c  oncentrations  s  uggesting  t  hat  a  pproximations  us  ed  i  n 
computing  t  he  w  ater  f  raction  w  ere  also  a  f  actor  i  n  t  hat  w  ork.  Most  of  t  he  i  mprovement  i  s 
probably  due  to  the  increase  in  the  quality  of  the  optical  hardware  over  the  course  of  the  last  8 
years  as  w  ell  as  m  ore  a  ccurate  d  etermination  o  f  t  he  ab  sorption  1  ine  pa  rameters.  S  ome  of  t  he 
differences  could  al  so  b  e  f  acility  i  mprovements.  H  owever,  i  t  i  s  al  so  t  he  cas  e  t  hat  t  he  u  se  o  f 
multiple  paths  allows  for  averaging  over  spatial  inhomogenities  that  would  be  more  difficult  with 
two  lines  of  sight.  The  work  of  Liu,  et  al.  represents  a  more  recent  study  (2005)  of  the  test  rig  in 
research  cell  22  a  nd  i  t  shows  t  he  i  mprovements  i  n  opt  ical  ha  rdware  t  hat  ha  ve  be  en  r  ecently 
made.  It  shows  many  of  the  improvements  in  data  quality  that  can  be  observed  by  comparing  the 
current  s  tudy  t  o  t  he  one  b  y  U  pschulte,  e  t  a  1.  U  nfortunately,  no  c  omparisons  t  o  the  pr  edicted 
facility  vitiation  conditions  were  made  so  a  straightforward  comparison  cannot  be  made  but  it 
appears  t  hat  t  he  m  easurement  q  uality  i  s  m  ore  i  n  ag  reement  w  ith  the  r  esults  g  iven  he  re.  T  he 
extent  of  LOS  multiplexing  (14  LOS)  and  data  acquisition  is  new. 
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The  de  sign,  s  etup,  a  nd  measurement  a  ccuracy  a  nd  pr  ecision  of  a  1 6  LOS  t  unable  di  ode  1  aser 
absorption  system  has  been  demonstrated.  Also,  the  first  ever  measurements  of  static  pressure 
using  1  inewidth  br  oadening  i  n  a  s  cramjet  i  solator  ha  ve  be  en  r  eported.  T  he  a  ccuracy  of  t  he 
measurements  is  excellent  with  average  differences  from  the  pressure  taps  on  the  order  of  2%  for 
static  p  ressure  and  t  he  temperature  a  greement  with  C  FD  b  eing  of  t  he  or  der  of  2%  .  T  he 
comparison  of  water  mole  percent  is  more  complex  and  requires  that  approximations  be  made  in 
facility  computations.  N  evertheless,  excellent  agreement  is  once  again  seen  with  the  exception 
of  the  run  AF  (x  ~  1.5%)  where  the  disagreement  is  17%  but  is  still  within  the  facility  error  bars. 
The  hardware  and  software  developed  for  this  study  pave  the  way  for  measurements  of  shock 
structures  in  the  isolator.  Furthermore,  the  system  described  here  is  scalable  and  should  enable 
high  pr  ecision,  s  patial  a  nd  s  pectrally  m  ultiplexed  T  DLAS  m  easurements  f  or  s  upersonic  flow 
characterization,  c  ombustion  pr  oduct  a  nalysis,  a  nd  r  eal-time  mo  nitoring  o  f  th  e  f  low 
environment.  At  present  16  beam  paths  are  employed,  but  the  system  could  be  readily  extended 
to  128  be  am  paths  with  the  appropriate  optical  access.  The  importance  of  this  advancement  is 
more  evident  Reference  6  which  focuses  on  the  analysis  of  transient  behavior. 


4.3  Accuracy  and  Precision  of  TDLAS  Measurements  in  Hypersonic  Flows 


The  s  ignificant  T  DLAS  m  easurement  ach  ievements  n  oted  ab  ove  p  oint  to  a  g  eneral  n  eed  f  or 
understanding,  in  a  detailed  way,  the  accuracy  and  precision  to  which  TDLAS  measurements  can 
be  made  in  reacting  flows.  While  t  he  r  elative  ac  curacy  and  p  recision  m  ay  v  ary  from  o  ne 

Table  6  Near-infrared  (-1.39  pm)  water  transitions  used  in  study. 


Transition 

Label 

Transition 

Frequency 

(cm1) 

Lower  Energy 

Level 

(cm1) 

Line  Strength  at  296  K 
(cm2  molecule  1  cm'1) 

A 

7161.410 

224.838 

1.168E-20 

B 

7179.752 

1216.195 

2.438E-22 

C 

7181.156 

136.762 

1.488E-20 

D 

7185.597 

1045.058 

8.335E-22 

particular  application  to  another,  the  limits  of  both  will  be  dictated  by  general  considerations  of 
signal-to-noise  ratios  and  signal  processing  algorithms.  In  applications  associated  with  reacting 
flows,  the  largest  source  of  uncertainty  is  typically  the  signal-to-noise  ratio  (SNR).  The  SNR  in 
turn  depends  on  intrinsic  detector  noise,  A/D  bit  depth,  parasitic  e talons,  gross  beam  steering, 
beam  obscuration,  etc.  Limits  of  accuracy  and  precision  will  also  be  detennined  by  the  accuracy 
of  the  fundamental  spectroscopic  co  nstants  as  sociated  w  ith  the  target  species  and  the  specific 
transitions  a  ccessed  dur  ing  s  ignal  c  ollection.  S  ome  of  t  his  t  hinking  h  as  be  en  addressed  a  s 
researchers  have  sought  resonant  features  that  exhibit  temperature  sensitivity  over  the  anticipated 
measurement  temperature  range.  We  are  currently  examining  such  issues  in  a  detailed  manner 
that  i  s  i  ndependent  o  f  h  ardware  d  etails  us  ing  T  DLAS  da  ta  acquired  i  n  a  num  her  of  flowing 
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environments.  T  he  following  di  scussion  i  s  focused  on  s  ingle-scan  (  un- averaged)  w  ater 
absorption  spectra  acquired  in  a  time-multiplexed  manner  over  the  transitions  shown  in  Table  6. 
The  spectra  were  acquired  using  direct  absorption  and  individually  processed  to  detennine  the 
integrated  area  under  each  spectral  feature.  T emperature  was  then  derived  from  the  integrated 
areas  using  two  approaches.  In  the  first,  ratios  of  integrated  areas  were  calculated  and  equated  to 
the  corresponding  ratio  of  the  line  strengths  as  shown  here 


A  i  _  LS  i 

~1  ~  ~ls~2 


(8) 


The  line  strengths  themselves  are  known  functions  of  temperature  and  the  ratio  of  two  is  simply 
given  by 


LS,  =  LS,  (T0 )  -i(£l-E2)(i-^) 

ls2  ls2(t0 ) 


(9) 


Shot  Number 


Figure  12.  Single-shot  TDLAS-based  temperature  measurements  (T  in  K)  acquired  in  a  Mach 

1.8  flow  stream  at  the  entrance  to  the  isolator. 

where  Et  denotes  the  lower  energy  of  the  ith  transition  and  T0  a  reference  temperature  (296  K). 
Extracting  a  te  mperature  v  alue  from  th  e  lin  e  ratio  me  thod  is  s  imply  a  ma  tter  o  f  finding  th  e 
unique  value  for  T  that  satisfies  Eqn.  8  and  9.  (The  line  strengths  at  the  reference  temperature  are 
found  in  the  open  literature.)  The  second  approach  for  determining  the  temperature  is  based  on 
the  w  ell-known  B  oltzmann  pi  ot.  In  t  his  p  rescription,  t  he  1  ogarithm  o  f  th  e  in  tegrated  areas 
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normalized  by  their  corresponding  lines  strengths  at  the  reference  temperature  are  graphed  as  a 
function  of  the  corresponding  values  of  the  lower-state  energy  level  for  the  transitions.  A  linear, 
least-squares  fit  to  the  plotted  data  yields  both  the  temperature  and  density  as  seen  from  Eqn.  A5 
in  t  he  Appendix.  T  he  de  rived  de  nsity  (N)  a  ssociated  w  ith  t  he  f  irst  a pproach  i  s  f  ound  di  rectly 
from  Aj  =  LSi(T)  N  L  using  the  temperature  as  found  from  Eqn.  8  and  Eqn.  9.  The  anticipated 
numerical  u  ncertainties  associated  w  ith  t  he  d  erived  t  emperatures  an  d  d  ensities  as  sociated  w  ith 
both  a  pproaches  a  re  w  ritten  e  xplicitly  i  n  t  he  Appendix.  B  oth  a  pproaches  pr  esume  1  ocal 
thermodynamic  equilibrium  of  the  target  species  over  the  time  duration  of  the  acquisition  of  a 
single  spectrum.  This  assumption  is  sound  under  the  conditions  that  the  data  presented  here  were 


Shot  Number 


Figure  13.  Single-shot  TDLAS-based  temperature  measurements  (T  in  K);  same  conditions  as 

Figure  12 

acquired. 

The  accu  racy  o  f  TDLAS-derived  q  uantities  s  uch  as  t  emperature  ca  n  b  e  as  certained  v  ia 
measurements  m  ade  i  n  s  teady-state  e  nvironments  w  ith  know  n  c  ontrolled  va  lues  or  i  n  s  imilar 
environments  i  n  w  hich  t  he  de  sired  pa  rameter  i  s  measured  us  ing  a  s  econd  know  n  m  ethod.  W  e 
have  acquired  water-based  TDLAS  measurements  (using  the  transitions  in  Table  6)  above  a  flat- 
flame  burner  along  with  corrected  thennocouple  measurements  in  the  same  methane/air  flames  at 
a  fixed  height  above  the  burner  surface.  (This  work  will  be  discussed  in  detail  in  a  future  paper.) 
For  the  particular  case  of  an  equivalence  ratio  of  0.9,  the  TDLAS  measurements  yielded  a  path- 
integrated  t  emperature  a  cross  t  he  bu  rner  of  182  5  K  us  ing  Boltzmann-plot  d  ata  an  alysis.  A 
sequence  oft  hermocouple  m  easurements  m  ade  with  t  hree  d  ifferent  b  ead  d  iameters  al  ong  t  he 
same  path  yielded  a  p  ath-integrated  temperature  value  of  1785  K.  These  two  values  agree  with 
each  o  ther  at  t  he  2  %  1  evel.  S  imilar  accu  racy  assessments  ap  pear  i  n  t  he  o  pen  1  iterature.  F  or 
example,  i  n  pr  eparation  f  or  w  ind- tunnel  t  esting  o  f  a  m  odel  s  cramjet  en  gine,  r  esearchers 
compared  water-based  TDLAS  measurements  of  temperature  with  those  of  a  thennocouple  in  a 
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premixed  flame.5  Over  the  temperature  range  o  f  1000  K  to  1600  K  t hey  found  t he  T  DLAS- 
derived  t  emperature  t  o  be  a  few  p  ercent  greater  t  han  t  he  t  hermocouple  t  emperature  w  ith  t  he 
difference  increasing  with  flame  temperature. 

For  me  asurements  ma  de  in  r  eacting  flows  -  particularly  t  hose  s  et  u  p  t  o  t  est  h  ardware  -  it  is 
frequently  not  pos  sible  t  o  obt  ain  a  s  econdary  temperature  m  easurement.  Indeed,  opt  ical 
diagnostics  are  typically  employed  for  temperature  measurements  in  environments  where  other 
thermometry  approaches  are  not  possible  or  practical.  It  is  prudent  therefore  to  examine  TDLAS 
data  reduction  strategies  with  regards  to  anticipated  numerical  scatter  given  finite  signal-to-noise 
ratios  in  the  context  of  apparent  measurement  scatter  (i.e.  precision).  Such  close  examination  can 
lead  to  conclusions  that  may  be  data-set  dependent  and  care  must  be  used  in  transferring  such 
considerations  between  data  sets. 

Using  transitions  A,  B  and  D  of  Table  6,  the  temperature  of  ambient  laboratory  air  was  measured 
following  the  Boltzmann  plot  and  two-line  approaches  for  extracting  temperature.  The  data  set 
included  1000  s  ingle  scans  of  the  employed  diode  lasers  acquired  at  1  k  Hz.  The  three  derived 
average  temperatures  agreed  to  within  6% .  ( We  note  that  the  transitions  of  Table  6  are  b  etter 
suited  f  or  elevated  t  emperature  e  nvironments  t  han  f  or  r  oom  t  emperature  e  nvironments.)  T  he 
Boltzmann  pi  ot  a  verage  t  emperature  o  f  293  K  a  grees  w  ith  t  he  t  wo-line  t  emperature  (288  K  ) 
using  the  B/A  line  ratio  at  the  1 .7%  level.  These  three  approaches  yielded  similar  average  values 
for  the  SNR-driven  average  uncertainty  of  ~33  K.  However,  the  observed  standard  deviations  of 
the  d  erived  t  emperatures  w  ere  s  mailer,  ~  2 1  K  .  T  he  1  arger  S  NR-driven  u  ncertainty  v  alue  i  s 
certainly  due  the  relatively  small  signal  associated  with  transitions  B  and  D  at  room  temperature. 

Table  7.  Extracted  temperatures  with  associated  uncertainties  using  Boltzmann  analysis 

and  area  ratios  for  transitions  of  Table  1. 


Boltzmann 

Plot  for 
A,B,C,D 

Boltzmann 

Plot  for 
A,B,C 

B/C 

B/A 

B/D 

A/C 

D/C 

D/A 

avg.  T  (K) 

488 

425 

427 

437 

300 

371 

479 

503 

std  expt.  (K) 

72 

49 

35 

38 

104 

122 

58 

68 

avg.  uncert.  (K) 

70 

84 

16 

30 

215 

198 

73 
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As  an  ex  ample  o  f  a  m  easurement  i  n  a  r  eacting  flow,  w  e  ex  amine  a  s  equence  o  f  s  ingle-scan 
(single-shot)  temperature  acquired  at  the  entrance  to  the  isolator  in  a  direct-connect  hypersonic 
test  facility.  The  measurements  were  made  using  a  single  line-of-sight  across  the  mid -plane  of 
the  circular  hardware  just  downstream  of  the  Mach  1.8  facility  nozzle.  The  incoming  air  stream 
consisted  of  heated  (non-vitiated)  air  with  an  average  water  molecule  number  density  of  6  x  10 16 
cm'3.  The  data  was  collected  at  1  kHz  using  the  transitions  noted  in  T able  6  and  processed  to 
yield  values  of  the  frequency- integrated  absorbance  for  each  shot.  Temperatures  were  derived 
from  these  integrated  areas  using  the  area  ratio  approach  (Eqn.  1  and  2)  and  the  Boltzmann  plot 
approach  (Eqn.  A  5).  The  f  our  t  ransitions  pr  ovide  f  or  s  ix  uni  que  a  rea  r  atios.  T  he  de  rived 
temperatures  for  an  arbitrarily  chosen  sequence  of  40  shots  are  shown  in  Fig.  12  and  13  for  the 
Boltzmann  plot  approach  and  five  of  the  ratio-based  approaches.  As  seen  in  Fig.  1 2,  the  results 
for  the  Boltzmann  pi  ot  approach  a  nd  t  he  r  atios  of  transition  c  ombinations  B/A  and  B/C(see 
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Table  6  for  transition  labeling  convention)  follow  each  other  quite  closely.  For  example  they 
each  exhibit  local  maxima  near  shot  numbers  23,  37,  and  52  and  local  minima  near  32  and  47. 
The  ap  parent  o  scillation  i  n  t  emperature  p  ersists  f  or  al  1 1  imes  r  ecorded  and  p  robably  r  eflects 
facility  air  supply  characteristics  such  as  compressor  surge.  Temperatures  derived  from  area  ratio 
combinations  involving  transition  D  (B/D,  D/C,  and  D/A)  exhibit  large  fluctuations  (see  Fig.  13) 
and  don’t  reproduce  the  periodic  structure  evident  in  Fig.  12.  Transition  D  exhibited  the  lowest 
signal-to-noise  r  atio  obs  erved  for  t  his  da  ta  s  et  and  i  s  t  he  1  ikely  r  eason  be  hind  t  he  excessive 
scatter  seen  in  Fig.  13. 

Using  a  full  second  worth  of  data  (1000  shots),  the  average  temperature  was  calculated  for  each 
approach  along  w  ith  t  he  s  tandard  de  viation  of  t  he  m  easurements.  T  hese  values  a  re  s  hown  i  n 
Table  7.  Similarly,  the  table  includes  the  average  of  the  uncertainties  of  the  derived  temperature 
values  based  on  m  easured  signal-to-noise  ratios  and  the  procedures  detailed  in  the  Appendix. 
Since  the  signal-to-noise  ratio  can  include  systematic  error  sources  we  do  not  add  the  individual 
uncertainties  in  quadrature  but  rather  report  the  simple  average.  Note  that  the  Boltzmann  plot 
derived  t  emperature  us  ing  all  f  our  t  ransitions  e  xhibits  1  arger  obs  erved  scatter  t  han  t  hat  of  t  he 
same  p  lot  with  tr  ansition  D  e  liminated.  T  he  te  mperature  d  erived  f  rom  th  e  th  ree-transition 
Boltzmann  analysis  is  plotted  in  Fig.  12. 

As  a  nticipated  from  F  igure  12  and  F  igure  13,  t  he  o  bserved  ex  perimental  s  catter  (standard 
deviation)  as  well  as  the  anticipated  numerical  scatter  are  least  for  the  three-transition  Boltzmann 
approach  and  the  ratio  combinations  of  B/C  and  B/A.  The  mean  temperatures  derived  via  these 
three  approaches  agree  within  2.8%.  The  aforementioned  oscillation  of  the  temperature  accounts 
for  the  measured  scatter  being  larger  than  the  anticipated  numerical  scatter.  The  two  area  ratio 
approaches  exhibit  1  ess  m  easured  s  catter  as  w  ell  as  1  ess  anticipated  numerical  u  ncertainty. 
Likewise,  as  anticipated  from  Fig.  13,  temperatures  derived  from  area  ratios  involving  transition 
D  exhibit  widely  varying  mean  values  along  with  large  measured  scatter  and  large  anticipated 
scatter. 

The  primary  quantities  derived  from  aTDLAS  measurement  are  the  temperature  and  number 
density  of  the  target  species.  Numerous  other  quantities  can  be  derived  as  well;  for  example,  the 
pressure  w  hich  i  s  de  rivable  f  rom  t  he  1  ine  w  idth  of  a  n  a  bsorption  t  ransition.  T  he  r  elevant 
expressions  a  re  s  hown  i  n  A  ppendix  Eqn.A12  -  A17w  here  u  se  h  as  b  een  m  ade  o  f  empirical 
relationships  between  the  full  width  at  half  maximum  of  the  transition,  the  temperature,  and  the 
pressure.  The  uncertainty  in  the  derived  pressure  is  then  dependent  upon  the  uncertainty  in  the 
derived  temperature  and  measured  line  width.  We  take,  as  an  example,  the  measured  temperature 
of  Table  7  (430  K)  and  the  Voigt  line  width  of  the  7185.6  cm'1  transition  (0.067  cm'1)  along  with 
constant  values  of  T o=296  K,  n  =  0.65,  a  nd  yo  =  0.041  c  nT  /atm  taken  from  the  literature  to 
find  a  measured  pressure  of  0.9  a  tm.  Assuming  no  unc  ertainty  in  the  empirical  constants,  the 
uncertainty  in  the  derived  value  of  the  pressure  is  0.07  atm  for  a  10%  temperature  uncertainty 
and  a  5%  line  width  uncertainty.  That  is,  the  combined  effects  of  temperature  and  line  width 
uncertainty  lead  to  an  8%  pressure  uncertainty. 

For  large  data  sets  with  thousands  of  single-shot  spectra  to  be  analyzed  it  is  helpful  to  develop 
routines  that  afford  some  means  of  rejecting  compromised  spectra.  Degradation  can  occur  for 
multiple  reasons  such  as  gross  beam  steering  or  obscuration  or  transient  electronic  noise.  One 
such  approach  that  we  have  found  valuable  is  close  examination  of  the  fitting  associated  with  the 
Boltzmann  plot.  The  scatter  about  the  linear  fit  can  be  quantified  using  the  standard  rules  for 
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calculating  %  .  For  a  single  data  set  at  constant  flow  condition  we  calculate  the  individual  values 
of  %  ,  average  them,  and  then  reject  outliers  that  show  large  deviation  from  this  average. 

While  i  ssues  as  sociated  w  ith  t  he  accu  racy,  p  recision,  an  d  s  catter  as  sociated  w  ith  f  low  f  ield 
parameters  d  erived  f  rom  T  DLAS  m  easurements  h  ave  b  een  d  iscussed  h  ere  u  sing  s  pecific  d  ata 
sets,  the  methodology  is  quite  general  and  can  be  extended  to  other  data  sets.  W e  recommend 
using  a  combination  of  Boltzmann  plot  analysis  combined  with  integrated  area  ratio  analysis  to 
extract  t  he  1  ocal  p  ath-integrated  te  mperature.  For  a  s  et  o  f  tr  ansitions  it  is  b  est  to  e  xclude  a 
transition  or  transitions  that  exhibit  poor  signal-to-noise  ratios  as  compared  to  other  transitions  in 
the  s  et.  A  s  i  ndicated  i  n  T  able  7  ,  cal  culation  of  t  he  a  nticipated  num  erical  s  catter  pr  ovides 
guidance  in  selecting  line  pairs  for  careful  analysis. 


5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  e  fforts  de  scribed  h  ere  1  ays  de  tailed  groundwork  f  or  t  he  us  e  of  T  DLAS  i  n  s  tudying  t  he 
dynamic  be  havior  of  scramjet  is  olators.  A  n  e  xample  o  f  d  etailed  me  asurements  o  f  shock  t  rain 
dynamics  was  presented.  Such  work  provides  invaluable  information  for  anchoring  CFD  codes 
and  related  fluidic  models.  As  yet,  there  is  no  c  omprehensive  theory  of  all  aspects  of  isolator 
perfonnance  and  behavior  preventing  detailed  prediction.  For  example,  it  is  not  now  possible  to 
accurately  p  redict  i  solator  m  argin  g  iven  v  alues  for  co  mbustor  h  eat  r  elease  an  d  p  ressure  fields. 
TDLAS  provides  a  tool  to  address  the  lack  of  fundamental  infonnation  which  will  lead  to  desired 
comprehensive  understanding. 

Future  work  should  strive  to  mature  the  TDLAS  technique  in  three  areas:  1)  benchmarking  of 
TDLAS  derived  physical  parameters  such  as  temperature  with  those  obtained  from  other  means, 
2)  broadening  and  refinement  of  the  basic  spectroscopic  database  for  species  such  as  water  and 
molecular  ox  ygen  to  i  mprove  t  he  a  ccuracy  of  a  vailable  s  pectroscopic  c  onstants,  a  nd  3) 
development  of  advanced  numerical  algorithms  for  extracting  spatiotemporal  infonnation  from 
TDLAS  m  easurements  m  ade  a  long  t  wo  or  m  ore  1  ines  of  s  ight.  Area  1  serves  t  o  f  acilitate  t  he 
future  us  e  o  f  T  DLAS  b  y  researchers  w  ho  do  n  ot  t  raditionally  us  e  s  pectroscopic  t  echniques. 
Benchmarking  TDLAS  measurements  against  other  sensor/metrology  standards  gives  such  users 
a  s  ense  o  f  ex  pected  p  erformance.  The  a  im  of  w  ork  i  n  a  rea  2  i  s  t  o  i  mprove  t  he  a  ccuracy  o  f 
parameters  derived  from  TDLAS  measurements.  The  ultimate  accuracy  cannot  be  better  than  the 
spectroscopic  parameters  used  in  the  TDLAS  data  processing  routines.  For  example  it  is  highly 
desirable  t  o  ha  ve  t  ables  of  a  bsorption  cross  s  ections  w  ith  i  ntrinsic  a  ccuracies  b  etter  t  han  1  % . 
Efforts  in  area  3  will  expand  the  utility  of  the  technique  making  comparison  to  modeling  efforts 
much  more  sophisticated  and  hence  valuable. 
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APPENDIX 


The  following  indicates  the  manner  in  which  the  numerical  uncertainties  in  the  derived  values  of 
temperature,  t  arget  s  pecies  num  her  de  nsity  and  pr  essure  di  scussed  above  a  re  found  us  ing 
standard  statistical  analysis  for  both  the  line  ratio  approach  and  the  Boltzmann  plot  approach. 

The  unc  ertainties  i  n  t  he  m  easured  a  bsorption  s  ignals  (A)  de  noted  by  a  A  are  d  etermined  f  rom 
signal-to-noise  conditions  as  noted  above.  For  any  function /that  is  dependent  upon  A  o r  a  set  of 
AjS,  the  uncertainty  in/is  detennined  from  the  following 


af2  = 


Al) 


where  z,  denotes  A,  or  a  function  of  A,.  Note,  we  are  presuming  independence  of  the  variables  z;-. 
Herein,  the  uncertainty  (cr)  is  expressed  as  a  standard  deviation.  The  95%  confidence  interval  can 
be  found  from  a  through  multiplication  by  1.96. 


We  first  consider  the  uncertainty  in  the  temperature  derived  from  a  two-line  ratio.  Defining  the 
ratio  of  the  integrated  absorbances  as  R  =  A1/A2,  the  fractional  temperature  uncertainty  can  be 
written  from  Eqn.  1  and  Eqn.  2  as 


aT 

T 


A  E 


gR 


fe/A E  Ln 


LS2(To ) 


)/lS1(To) 


+  1/to 


A2) 


where  AE  =  Ej  -  E2.  Similarly,  the  uncertainty  in  the  extracted  number  density  can  be  expressed 
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A  Boltzmann  analysis  of  absorption  spectra  stems  from  the  following  expression  for  each  (ith) 
transition  in  th  e  s  pectrum  (  we  o  mit  f  actors  a  ssociated  w  ith  s  timulated  e  mission  th  at  a  re 
insignificant  in  this  context) 
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where  To  denotes  a  reference  temperature  (usually  296  K  ).  A  Boltzmann  plot  is  generated  by 
graphing  the  left  hand  side  of  Eqn.  A5  as  a  function  of  Et  for  two  or  more  transitions.  A  linear  fit 
to  the  plotted  points  yields  the  temperature  ( T)  via  the  slope  and  the  number  density  (N)  of  the 
target  species  via  the  y  intercept.  To  extract  a  measure  of  the  uncertainty  in  the  derived  values  for 
T  and  N,  we  write  the  linear  fit  to  the  Boltzmann  plot  in  the  form  y  =  a  +  b  x.  Uncertainties  in  a 
and  b  are  derivable  from  the  uncertainties  in  the  values  of  y  that  in  turn  reflect  the  uncertainties 
in  the  set  of  At.  A  solution  for  a  and  b  is  found  using  well-established  routines  that  we  reproduce 
here  following  the  notation  of  Press  et  al. 1  Each  value  of  v*  has  an  associated  uncertainty  noted  as 
Of- 
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The  fitted  values  of  a  and  b  along  with  their  uncertainties  are  then  given  by 
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Following  a  Boltzmann  analysis,  the  unc  ertainty  in  the  derived  value  for  the  temperature  and 
target  species  number  density  are  given  from  the  above  by 
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1  Press,  W.  H.,  Teukolsky,  S.  A.,  Vetterling,  W.  T.,  Flannery,  B.  P.,  Numerical  Recipes  in  C++,  Chap.  15, 
(Cambridge  University  Press,  Cambridge,  England,  2002). 
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TDLAS  measurements  can  be  used  to  determine  the  local  pressure  through  the  line  widths  (full 
width  at  half  maximum)  of  the  recorded  spectral  transitions.  Under  the  physical  environments 
associated  with  hydrocarbon-based  reacting  flows,  the  spectral  line  shapes  are  well  characterized 
by  the  Voigt  spectral  profde.  This  profde  includes  line  broadening  due  to  temperature  (Doppler 
broadening)  as  well  as  pressure  (collision  broadening).  The  collision  broadening  contribution  to 
the  line  width  can  be  separated  from  the  total  line  width  and  Doppler  contribution  as  follows 
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where  fi  denotes  c  onstants  a  ssociated  w  ith  D  oppler  w  idth  expression  ( AvD  =  [1  T/:).  C  areful 
experimental  work  has  revealed  scaling  laws  that  empirically  show  the  relationship  between  the 
collision  width  and  the  pressure  that  are  often  written  in  the  following  fonn 
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where  n  and  y0  are  transition-specific  constants  and  T0  the  reference  temperature.  The  above  can 
be  combined  into  a  single  expression 
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The  statistical  uncertainty  in  the  pressure  can  then  be  found  from 
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